Dynamical properties of a 3-dimensional diluted Heisenberg model. 
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We study the magnetic excitation spectrum in 3D diluted ferromagnetic nearest neigbour sys- 
tems down to the percolation threshold. The disorder effects resulting from the dilution are handled 
accurately within Self-Consistent Local RPA approach. The calculations are performed using rel- 
atively large systems containing typically 20000 localised spins, a systematic average over many 
configurations of disorder is performed. We analyze in details the change in the magnon spectrum 
and magnon density of states as we increase the dilution. The zone of stability of the well defined 
magnon modes is shown to shrink drastically as we approach the percolation threshold. We also cal- 
culate the spin stiffness which appear to vanish at the percolation threshold exactly. A comparison 
with available data, based on a different theoretical approach, is also provided. We hope that this 
study will motivate new experimental studies based on inelastic Neutron Scattering measurements. 



Over the last few decades interest in disordered 
magnetic systems, like transition metal alloysi^— , di- 
luted magnetic semiconductors as Mn doped III-V 
compounds^"—, dP matcrialsi^"— and manganites^i"— 
has considerably increased. Most of these materials 
have attracted a great deal of attention from both the 
experimental as well as the theoretical point of view 
due to their potential in spintronics applications. In- 
elastic neutron scattering is the most direct experimen- 
tal method to study the magnetic excitations in these 
materials. However, theoretical studies require tools 
which can treat both disorder and thermal fluctuations 
in a reliable manner. Different theoretical methods have 
been performed/developed to calculate magnetic prop- 
erties and the transition temperatures. Among the ap- 
propriate tools one can first quote the classical Monte 
Carlo method which is essentially exact. However to our 
knowledge no studies on the magnetic excitation spec- 
trum in dilute 3d Heisenberg model has been performed, 
probably because of the cost both in memory and CPU 
time. Most of the existing studies within classical Monte 
Carlo method deal with the equilibrium properties of 
the Ising model^^"— and there are a few focusing on the 
3D diluted nearest neighbor Heisenberg model2^. There 
is an alternative approach based on finite temperature 
Green's functions known as Self-Consistent Local Ran- 
dom Phase approximation (SC-LRPA). It was shown to 
be reliable to handle dilution/disorder effects and ther- 
mal fluctuations^"—. In the SC-LRPA theory disor- 
der/dilution is treated exactly while thermal fluctuations 
are treated within random phase approximation. This 
method has several advantages, it is semi-analytical, fast, 
and allows to reach very large systems. The aim of the 
present manuscript is to provide a detailed study of the 
magnetic excitation spectrum in the dilute ferromagnetic 
nearest-neighbor Heisenberg Hamiltonian. 



formalism. The Hamiltonian describing Nimp interacting 
quantum/classical spins randomly distributed on a lat- 
tice of N sites is given by the diluted random Heisenberg 
model, 



(1) 



where the couplings Jij are assumed to be given and are 
restricted to nearest neighbor only (J). The sum ij runs 
over all sites and the random variable p; is 1 if the site 
is occupied by a spin otherwise it is 0. In the present 
work, the calculations are performed at T = K, the 
quantum/classical nature of the spin is irrelevant. Thus 
wc assume classical spins with | Si | = 1. 

We introduce the retarded Green's func- 
tions G^/c.)=/_+^ G^/t)e-*dt where G'r^{t)= - 
i6{t){[S^ {t), S~ {0)]), which describe the transverse spin 
fluctuations and the subscript 'c' corresponds to the con- 
figuration of disorder. After performing the Tyablicov 
decoupling of the higher order Green's functions which 
appear in the equation of motion of G^ ■ (w)^, we obtain 



M-H^„)G^ = D 



(2) 



where H^y^jG^ and D are Ni^p x Ni^p matrices. 

The effective Hamiltonian matrix elements are 
(H^^^),, = -{S^} J,, + <5,,Ez(^z^)Jy and D„=2(^f )<5„-. 
For a given temperature and disorder configuration, 
the local magnetizations {S^) {i = 1,2, ....Nimp) have 
to be calculated self-consistently at each temperature. 
The local magnetization is evaluated from the Callen 
expressions^, which relates the local Green's function at 
site i to the local magnetization at this site, 

{s ^iKi + <^^?'+' + ii + s + $0$.^^+^ 
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The local effective magnon occupation number is given 



Let us now present a short summary of the SC-LRPA by $j — 
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FIG. 1: (Color online) Magnon density of states p{uj) as a 
function of the energy cj for different concentration of localised 
spins x.The energy axis (x axis) is in oj/J. 



Despite the matrix being non-Hermitian, the spec- 
trum is real and positive at each temperature in the 
ferromagnetic phase. A negative eigenvalue would in- 
dicate an instability of the ferromagnetic phase. H^^^ is 
non-Hermitian (real non-symmetric) but it has the prop- 
erty to be biorthogonal^i. Hence one has to define the 
right and left eigenvectors of H'^fj denoted by I'^^a''^) ^-^d 
'I'^'^) respectively, both associated to the same eigen- 
value ll)^. After inserting the L and R eigenvectors in 
eq.(2), the retarded Green's functions can be written as 
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The averaged Fourier transform is given by G'(q, w) = 

(lvl!;:;E»,e"'^''""'''^G^^('^))c: where the notation ( ), 

denotes the average over the disorder configurations. 
The dynamical spectral function is given by A(q, cj) = 
-^^ImG(q,w), where {{S^}} = j^E^iSf) is the 
total magnetization averaged over all spin sites. The 
spectral function A{q,uj) provides direct access to the 
magnetic excitation spectrum, it is also directly accessi- 
ble by Inelastic Neutron Scattering experiments. 

One can also calculate the magnon density of states 
which is given by p{uj) = J2iPii^)^ where pi{uj) = 
— 27^^ImGii(a;), is the local magnon DOS. Here we 
have calculated the magnetic excitation spectrum for a 
simple cubic lattice as a function of the concentration of 



the spins a; 



N 



The calculations arc performed at 



T=0 K, where the matrix Heff is real symmetric and 
hence L and R eigenvectors arc identical. Since the cou- 
plings used in this case are all ferromagnetic, there is no 
frustration and the system exhibits long range ferromag- 
netic order beyond the percolation threshold. 

In Fig[T]we have plotted the magnon density of states 
(DOS) as a function of the energy, for different concentra- 
tions of localised spins {x—1.0 corresponds to the clean 



case). As the dilution increases, we observe (i) a signifi- 
cant change in the shape of the magnon DOS, (ii) a re- 
duction of the bandwidth and (iii) an increase in the low 
energy DOS. As we approach the percolation threshold 
which is located at a;c=0.31— , we observe for the low- 
est concentration a sharp peak developing at a;=0. This 
feature is attributed to the formation of several isolated 
clusters which have their own zero energy modes. Addi- 
tionally we observe at higher energy, peak structures lo- 
cated at w= J and 2 J respectively. These peaks are due to 
the non-zero eigenmodes of these isolated clusters. Note 
that the DOS does not provide any information on the 
nature (extended/localized) of the magnon modes. This 
would require a more careful analysis and study, either 
using the calculation of the inverse participation ratio 
(IPR) or the typical magnon DOS^ which provide a di- 
rect access to the "mobility edge" separating the localized 
modes from the extended ones. The typical magnon DOS 
is given by ptypitj) = exp ((In Pi{uj))J , where ( )c rep- 
resents the average over disorder configurations and Pi{ui) 
is the local magnon DOS. We expect the mobility edge 
to become zero at the percolation threshold exactly, for 
a review see 
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FIG. 2: (Color online) Spectral function j4(q, a;) as a function 
of the energy uj in the (1,0,0) direction for different values 
of qx, with varying concentration of localised spins (a::).The 
energy axis (x axis) is in lj/ J 

FigUl shows the spectral function A(q, a;) as a func- 
tion of energy for different values of the momentum q in 
the (100) direction. We have shown the spectral func- 
tion for different concentrations of spins a;, with x~1.0 
corresponding to the clean case. The localised spins arc 
randomly distributed over a simple cubic lattice of size 
(L=28) and we have used periodic boundary conditions. 
For x=0.8, the system contains approximately 17000 
spins. Wc have examined the finite-size effect by compar- 
ing with the results for L=20 and L=24, and concluded 
that L=28 is large enough for the following discussions. 
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The average over disorder was done for 50 configurations. 
We have also checked for the necessary number of ran- 
dom configurations and found that 20 configurations are 
sufficient for L=28 system. The number of disorder con- 
figurations should be sufficient in order to extract the 
correct spin stiffness. We can clearly see from the fig- 
ure that well defined excitations exist only for small val- 
ues of q. As the momentum increases the peak becomes 
broader and develops a tail extending towards higher en- 
ergies. As we get closer to the percolation threshold, the 
magnetic excitation peaks become strongly asymmetric. 
Similar asymmetric peaks were also observed in the mag- 
netic excitation spectrum for the diluted antiferromagnet 
MnxZni_a;F2^. This increase in asymmetry with the 
momentum corresponds to a crossover from propagating 
spin waves to localized excitations (fractons) '^^'^" . 
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FIG. 3: (Color online) Spectral function ^(q, oj) in the {q, u) 
plane for different concentration of localised spins (x). (The 
size of the simple cubic lattice is L=32) . 

In the next figure, FiglU we have plotted the spectral 
function in the {q, lu) plane, corresponding to three dif- 
ferent concentrations of spins. The spectral function is 
shown for the entire Brillouin zone. Well defined exci- 
tations can be observed in the dilute case only in a re- 
stricted region of the Brillouin zone close to the F point 
{q={000)). Similar results were observed for the exci- 
tation spectrum of the diluted magnetic semiconductor 
Gai_2;Mn2;A9^. From the figures we conclude that well 
defined excitations cease to exist as we go to higher di- 
lutions and the spectrum shows broadening close to the 
percolation threshold due to the formation of localized 
modes. 

The spin wave energy uj{q) for different concentrations 
of spins as a function of is shown in FigH) Actually 
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FIG. 4: (Color online) Magnon energy uj{q) (in uj/J) as a 
function of for different concentration of localised spins (x). 
The size of the simple cubic lattice varies from L=20,24,28 and 
32. 




FIG. 5: (Color online) Spin stiffness D (in D/J) as a func- 
tion of the concentration of the localised spins (x). The dia- 
monds represent the data from ref.— .The squares are calcu- 
lated within SC-LRPA. The dashed line is the Virtual Crystal 
Approximation value, (sc is the percolation threshold). 

the energy of the first peak in ^(q, w) in the (100) direc- 
tion is shown as a function of q^. The system size varies 
from L=20 to L=32. As can be seen from the figure for 
small momentum, uj{q) « D{x)q^. The slope of these 
curves gives the spin stiffness D{x) for different x. We 
observe a strong decrease of the slope as we approach 
the percolation threshold. Note that, the perfect linear 
behaviour of w(g) as a function of vindicates the fact 
that the average over disorder configurations was obvi- 
ously sufficient. In FigOwe have plotted the spin stiffness 
D{x) as a function of the concentration of the localised 
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spins X. We find that the spin stiffness is almost hnear 
from x = 1.0 down to x = 0.5 and only forms a concave 
toe close to the percolation. We have also shown the 
spin stiffness obtained by Kirkpatrick^.The author has 
used site percolation statistics on random resistor net- 
works, as a natural generalization of lattice models, and 
obtained a relation between the site percolation proba- 
bility P^*'(a;), the conductance G{x) and the spin stiff- 
ness coefficient D{x). We find that our results are in 
very good agreement with those of refi^. Note that, the 
method we have used here is entirely different since the 
spin-stiffness is extracted directly from the curvature of 
the magnons excitations. We have also plotted the Vir- 
tual Crystal Approximation (VGA) value which is given 
by D^'~^^{x) = xD{x = 1). As we decrease the concen- 
tration of localized spins the difference between the VGA 
and the SG-LRPA value increases significantly. Below 
X = 0.5 the D^'~'^{x) drastically overestimates the cal- 
culated values. Note that our spin stiffness curve is very 
different from that obtained in the case of optimally an- 
nealed samples of the diluted magnetic semiconductor 
Gai-ajMuj; As^. This difference of behavior is due to the 
fact that in the case of Gai_2^Mna;As the concentrations 
were much smaller and that the couplings were rather 
extended and spin concentration dependant. 

In conclusion, we have calculated the magnetic excita- 
tion spectrum for a dilute and disordered ferromagnetic 



system by using an approach based on the SG-LRPA, 
which allows us to treat disorder and thermal fluctuations 
in a reliable manner. We see that a proper treatment of 
the disorder leads to unusual excitation spectrum com- 
pared to what is usually observed in non dilute systems. 
We have also calculated the spin stiffness, without adjust- 
ing parameters, as a function of the spin concentration 
and observed that our results are in very good agreement 
with previous studies for a similar model. Earlier the SG- 
LRPA was proved to be reliable for the determination of 
the Gurie temperature^. Here we have shown its accu- 
racy and efficiency for the calculation of the magnetic 
excitation spectrum. We hope that this study will stim- 
ulate new experimental measurements like the dynamical 
spectral function and spin stiffness as a function of the 
localised spin concentration. 
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